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We develop a nonperturbative method to calculate the electron propagator in low-dimensional 
and strongly correlated electron systems. The method builds on our earlier work using a Hubbard- 
Stratonovich transformation to map the tunneling problem to the x-ray edge problem, which ac- 
counts for the infrared catastrophe caused by the sudden introduction of a new electron into a 
conductor during a tunneling event. Here we use a cumulant expansion to include fluctuations 
about this x-ray edge limit. We find that the dominant effect of electron-electron interaction at 
low energies is to correct the noninteracting Green's function by a factor e~ s ^ h , where S can be 
interpreted as the Euclidean action for a density field describing the time-dependent charge distri- 
bution of the newly added electron. Initially localized, this charge distribution spreads in time as 
the electron is accommodated by the host conductor, and during this relaxation process action is 
accumulated according to classical electrostatics with a screened interaction. The theory applies 
to lattice or continuum models of any dimensionality, with or without translational invariance. In 
one dimension the method correctly predicts a power-law density of states (DOS) for electrons 
with short-range interaction and no disorder, and when applied to the solvable Tomonaga-Luttinger 
model, the exact DOS is obtained. 

PACS numbers: 71.10.-w, 73.43.Jn 



I. INTRODUCTION 

In the conventional many-body theory treatment of 
tunneling, the low-temperature tunneling current be- 
tween an ordinary metal and a strongly correlated elec- 
tron system is controlled by the single-particle density of 
states (DOS) of the correlated system. Tunneling exper- 
iments are therefore often used as a probe of the DOS. 
In a wide variety of low dimensional and strongly corre- 
lated electron systems, including all ID metals, the 2D 
diffusive metal, the 2D Hall fluid, and the edge of the 
sharply confined Hall fluid, the DOS exhibits anomalies 
such as cusps, algebraic suppressions, and pseudogaps 
at the Fermi energy. We have recently proposed that 
the origin of these anomalies is the infrared catastrophe 
caused by the response of the host electron gas to the 
sudden introduction of a new particle that occurs during 
a tunneling event i 

The infrared catastrophe is a singular screening re- 
sponse of a degenerate Fermi gas to a localized poten- 
tial applied abruptly in time, which is known to be re- 
sponsible for the anomalous x-ray optical and photoemis- 
sion edge spectra of metals the Anderson orthogonal- 
ity catastrophe,^ and the Kondo effectual Imagine an 
electron in a tunneling process being replaced by a nega- 
tively charged, distinguishable particle with mass M, In 
the M — >oo limit, the time-dependent potential produced 
by the particle being added to the origin at time r = 
and removed at a later time To is& 

<Mr,r) =f/(r)e(r)e(r -r) ) (1) 

where U is the two-particle interaction. xr is identical, 
up to a sign, to the abruptly turned-on hole potential 
of the x-ray edge problem, so an infrared catastrophe 



would be expected. Tunneling of a real, finite-mass elec- 
tron is different because it recoils, softening the potential 
produced. The real electron is also an indistinguishable 
fcrmion, unable to tunnel into the occupied states below 
the Fermi energy. 



In Ref. Q we introduced an exact functional-integral 
representation for the interacting propagator, and devel- 
oped a nonperturbative technique for calculating it by 
identifying a "dangerous" scalar field configuration of the 
form Q , and then treating this special field configuration 
by using methods developed for the x-ray edge problem. 
All other field configurations were ignored, thereby re- 
ducing the tunneling problem to an x-ray edge problem. 
Nonetheless, qualitatively correct results were obtained 
for the ID electron gas and the 2D Hall fluid using this 
approach. In this paper we attempt to go beyond this 
so-called x-ray edge limit by including fluctuations about 
cj> XI . We find that by including fluctuations through the 
use of a simple functional cumulant expansion, a qualita- 
tively correct DOS is obtained for electrons with short- 
range interaction and no disorder in one, two, and three 
dimensions. We also show that and when applied to 
the solvable Tomonaga-Luttinger model, the low-energy 
fixed-point Hamiltonian for most ID metals, the exact 
DOS is obtained. A preliminary account of this work 
has appeared elsewhere^ 
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II. GENERAL FORMALISM AND CUMULANT 
EXPANSION METHOD 

We assume a grand-canonical Hamiltonian H = Ho + 
V, where 



III 

2m 



+ u(r) - fx 



Mr) (2) 



is the Hartree mean-field Hamiltonian, and 



V = i / d D r dV Sn(r) U(r - r') <Jn(r'). 



(3) 



Here II = p + |A, and u(r) is any single-particle po- 
tential, which may include a periodic lattice potential or 
disorder or both, and which also includes the Hartree 
interaction with the self-consistent density no(r). The 
interaction term is written in terms of the density fluc- 
tuation 



5n(r) = ^^t(r)<i/v(r) - ra (r). 



(4) 



Following Ref. l|, we use a Hubbard-Stratonovich trans- 
formation to write the exact Euclidean propagator 

G(r f cr f ,r i cr i ,r ) = -(Ttp a! (r f , T O )-0 CTi (ri, 0))^, (5) 

in the form 



G(rfCT f ,ricri,ro) = N ■ 



where 



jD<j>e-y 



(, r f cr f ,ricri,ro| 



(6) 



g{4>) = -(T^ ( (r f) r )^ i (r i ,0)e i /^ r ^ <5 "^)> (7) 

is a noninteracting correlation function in the presence 
of a purely imaginary scalar potential i(f>(r, r), and TV = 

(Texp(— J drV))^ 1 is a constant, independent of To. 
Eq. © is an exact expression for the interacting Green's 
function. 

The region of function space that contributes to the 
functional integrals in © is controlled by the width of the 
Gaussian, which in the small U limit becomes strongly 
localized around = 0. By expanding J7J in powers 
of <j) and doing the functional integrals term by term, 
one simply recovers the standard perturbative expansion 
for G(rfCTf , r;CTi, r ) in powers of U. Therefore, it will 
be necessary to go beyond a perturbative expansion for 
<7(rf<Tf, rj<Ti, To\<p). We evaluate ^(rfcrf, r;CTi, To|(/>) approx- 
imately, using a second-order functional cumulant expan- 
sion. Such an expansion amounts to a resummation of 
the most divergent terms in the perturbation series when 
(j) = cj) XI and the infrared catastrophe occurs. Indeed, one 
can view our resulting expression for g(rtat, ricTi, tq\4>) as 
a functional generalization of Mahan's "perturbative" re- 
sult for a similar correlation function. 2 Furthermore, for 
field configurations far from xr , the cumulant expansion 



will yield a result that is, by construction, exact through 
second order in U. After carrying this out we obtain 



g{4>) « Go(r f cTf,ricri,To) 



x e 

where 
and 

G 2 (rr,rV) = 



/ Ci(rr)0(r,r)+J C 2 (rr.rV) </>(r,r) 0(r' ,r') 



, (8) 



Gi (it) 



9i (rr) 



Go(rf<7f,ri<7i,To) 
g 2 (rr,rV) g x (it) g x (rV) 



(9) 



G (r f crf,r i cr i ,T ) 2 [G (r f (7f , ncri, r )] 2 ' 

(10) 



Here 

5n(riri, ■ • • , r„r„ 



y(TV' CTf (r f ,To)?A (Ti (ri,0)(5n(riri) ••• Sn(r n T n )) Q 

(11) 



is the coefficient of 4> n appearing in the perturbative ex- 
pansion of 10, as in 

5>(r f c7f,ricri,To|0) = G (r f cr f , r^, r ) 

oo - 

+ / 5r»(rm,--' ,r„T„)0(ri,ri)---0(r n ,r„). (12) 

71=1 ^ 

The cumulants Gi and G2 in terms of Go are 

„, \ G (r f cr f ,rcr, r -r)Go(rcr,riCTi,r) 

Oi(rr) = -i > — 

^ Go(rf(Tf,ricri,ro) 

(13) 

and (suppressing spin for clarity) 

C2(rr,rV) = -— 

2Go(r f ,ri,T ) 2 

x |[G (r,r',r - r')G (r f , r b r ) 

- G (r, n, r)G (r f , r', r - r')][r ^ r', r ^ r'] } (14) 

The functional integral in © can now be done exactly, 
leading to 

G(rf<7f,ri<7i,T ) =A(t ) Go(Tta { ,Tiai,Ta) e" s(ro) , (15) 
where 

f L>0 e -f J* -Kf^-ac^ 



/V[det(l -2G 2 C/)]~*, 



(16) 



and 



5=1 J^drdr' fd D rd D r'p(r,T)U eS (rT,r'T')p(r',T'). 

(17) 
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Here 



p(r,r) = -iC7i(rT) (18) 
Go(r f crf,rcr,T - r) G (rcr, riffi, t) 



G (r f tTf,ri(Ti,To) 



(19) 



and 



U eS (rr, r'r') = [ET^r - r')5(r - r') - 2G 2 (rr, r'r')] 1 

(20) 

is a screened interaction. 

Because spin-orbit coupling has been neglected in H 7 
the noninteracting Green's function is diagonal in spin, 
and 

. . Go(r f o'i,ro'i,ro -T)Go(ro-i,ricri,r) 

P\ V :T) = — - d CTiCTf (21) 

Go(rfO-i,ri«7i,T ) " " 

Eq. I|15|) is the principal result of this work. 



III. CHARGE SPREADING INTERPRETATION 

We interpret (fH))) as follows: S is the Euclidean 
action 10 for a time-dependent charge distribution p(r, r). 

I 



We shall show that p(r, t) acts like a charge density as- 
sociated with an electron being inserted at r; at t = 
and removed at r-f at to. This charge density interacts 
via an effective interaction U e g(rT, r'r') that accounts 
for the modification of the electron-electron interaction 
by dynamic screening^ Our result can therefore be re- 
garded as a variant of the intuitive but phenomenologi- 
cal "charge spreading" picture of Spivaki^ and of Levitov 
and Shytov^ However, here the dynamics of p(r, r) is 
completely determined by the mean-field Hamiltonian, 
and has the dynamics of essentially noninteracting elec- 
trons. 

First consider the integrated charge, 



Q{r) 



d D rp(r,r) 

J d D r G (r f ,r,To - T)G (r,ri,r) 
Go(r f ,ri,To) 



(22) 



Using exact eigenfunction expansions for the Green's 
functions we obtain 



$;(r f )$ a (r i )e-( e °-^ T « {[n F (e Q - M ) - l] 2 6(r - r)9(r) + n F (e a - p)[n F {e a - M ) - 1] [6(-r) + 6(r - r )]} 



0(7 



$>:(r f )* a ( ri ) ( 



-(e a -n)T 



n F (e a - fi) - 1] 



where n F is the Fermi distribution function and the <i> Q are the single-particle eigenfunctions of Hq. In the zero 
temperature limit 

K(*f)* a (*i)e~ lea ~" )n {( N * - i) 2 ©^ - t)Q(t) + N a (N a - 1) [9(-r) + 9(r - r )]} 

Q(r) = — ^ , (24) 

^<f>;(r f )$ Q (r i )e- (e °^ )T0 (^a-l) 

a 

i 



where N a is the ground-state occupation number of state 
a, which in the absence of ground state degeneracy takes 
the value of or 1. In this case reduces to 



Q(r) = O(r - r)6(r). 



(25) 



When the sum rule (|25|l holds, the net added charge, 
as described by p(r, r), is unity (in units of the electron 
charge) for times between and to, and zero otherwise. 
This behavior correctly mimics the action of the field 
operators in (JSJ. 

At short times, r < To, the charge density is approxi- 
mately 



which is localized around r = r;. As time evolves this 
distribution relaxes. Then as r approaches To the charge 
density again becomes localized around r = rf , 



p(r,r) Go(r f ,r,T - r). 



(27) 



A plot of p(r, t) for the ID electron gas is given in Fig.^ 
The dynamics of p(r, r) can be shown to be governed 
by the equation of motion 

[a T +iT (r)]p(r,T) = -<5(r f -r)S(r -r') + 6(t- Ti )6(t). 

This can be seen by noting that the noninteracting 
Green's function satisfies 



p(t,t) Go(r,r;,r) 



(26) 



[d T + H (r)} G (r, r', r, r') = -S(r - t')S(t - t'). (29) 
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FIG. 1: Charge density p{x, r) for the ID electron gas at two 
times, showing Friedel oscillations and gradual spreading. 

Then, using the definition (|19|l one can obtain (|28J) . 
Again, we stress that p(r, r) is describing the dynamics 
of noninteracting electrons, governed by the mean-field 
Hamiltonian Hq. 

Although p(r, t) has many properties that make it 
reasonable to interpret as the charge density associated 
with the added and subsequently removed electron in the 
Green's function JSJ), one should not take this interpre- 
tation too literally. For instance, the sum rule (|25|l only 
holds in the zero temperature limit and in the absence of 
ground state degeneracy. Also, as will be seen below in 
Sec. IIVCI p(r, t) may even be complex valued. 

IV. APPLICATIONS OF THE CUMULANT 
METHOD 



p to vary as £ . The interaction energy of such a charge 
distribution (assuming a short-range interaction) is 

E =\\ d ^[pW] 2 (30) 

which varies with time as t~ d . The action accumulated 
up to time To therefore scales as 

S ~ — nZTT if D > 2 ( 31 ) 

or 

S^UIutq if D = l. (32) 

The cases (|31[) and (|.':>2t arc dramatically different: In 
2D and 3D the action vanishes a long times, and the 
propagator (|T5j) is therefore not appreciably affected by 
interactions. The resulting DOS is energy independent 
at low energies, and the expected Fermi liquid behavior 
is recovered. In ID, however, the action diverges loga- 
rithmically, leading to an algebraic DOS. 

B. ID electron gas: Recovery of the Luttinger 
liquid phase 

The scaling argument of the previous section showed 
that the DOS in the ID electron gas with short-range 
interaction is algebraic, as expected. In this section we 
calculate the associated exponent. 

We proceed in two stages. Initially we keep only the 
first cumulant Ci, and then afterwards we discuss the 
effect of CV In ID it is possible to calculate the action 
(|17[) exactly in the long-time asymptotic limit at the first- 
cumulant level. Setting Xf = x x = 0, we have [see <|15[l ] 



In the following examples we will assume electrons with 
a short-range interaction U , no disorder, and no magnetic 
field. In Sec. IIV Al we show that our method correctly 
predicts a constant DOS near the Fermi energy in 2D 
and 3D, and in Sec. IIV _BI we obtain a power-law DOS 
in ID, in qualitative agreement with Luttinger liquid 
theory 14 ! 15 ! 16 ! 17 . 18 ! 19 ! 20 Finally, in Sec. ITVCl we use our 
method to calculate the DOS for the solvable Tomonaga- 
Luttinger model, obtaining the exact DOS exponent. 



A. 2D and 3D electron gas: Recovery of the Fermi 
liquid phase 

The sum rule (|25(l allows one to determine the energy 
dependence of the DOS in D dimensions, asymptotically 
in the low-energy limit, as follows: In the absence of 
disorder, a droplet of charge injected into a degenerate 
Fermi gas with velocity vp will relax to a size of order 
£ ~ v-pT after a time r. Approximating p(r, r) to be of 
uniform magnitude in a region of size £ and vanishing 
elsewhere, the sum rule then requires the magnitude or 



G(t ) = const x G (t ) e - 5(ro) . 



(33) 



Considering a local interaction of the form U(x — x') 
UqXS(x — x') the action is 



S ( T o) = ~7r~~ [ dr f dx [p{a 

1 JO J -oo 



(34) 



By linearizing the spectrum around the Fermi energy, the 
zero-temperature propagator at low energies is 



Gq{x,t) = — Im 

7T \X + WFT 

x sin k^x — v-pT cos k-px 



w(x 2 + v f t 2 ) 
The charge density (|19fl in this case is 

vpTo x sin kpx — vf(tq — r) cos kpx 



(35) 



p(x,r) 



7T X 2 + Vp(T — T ) 2 

x sin kpx — vpT cos k-px 

x 2 + V%T 2 



(36) 
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Fig. shows the charge density p(x, r) as it spreads in 
time from its initially localized position. 

By a lengthy but straightforward calculation it can be 
shown that 



S(t ) 



3U q X 



hi 



vftt 



(37) 



where a is a microscopic cutoff. This leads to a power-law 
decay of the interacting propagator as 

G(t) ~ ± (38) 



where 



a = — - 



3C/nA 



(39) 



8 VpTT 

is the propagator exponent. The DOS exponent 6, de- 
fined as 



N(e) = const x e s 



is given in this case by 



5 _ 3 U Q X 

8 VFTT 



(40) 



(41) 



The effects of the second cumulant C 2 are now straight- 
forward to understand: In addition to introducing a 
slowly varying prefactor A, whose only To dependence 
comes from the time dependence of the screening in (|20H , 
the second cumulant screens the bare interaction and 
does not prevent the logarithmic divergence of the ac- 
tion, but it does modify the DOS exponent. It is inter- 
esting, however, that in the large U limit, the effective 
interaction becomes independent of U, a clear indication 
of nonperturbative behavior. 

It is illustrative to compare (|41|l to the prediction of 
the perturbative x-ray edge limit of Ref. U where one 
neglects all field configurations in © except 4> XI . There 
we found (for this same short-range interaction model), 



U X 



vftt 



(42) 



to leading order in U. The DOS exponents l|41|) and i|42|) 

are in qualitative agreement, but the inclusion of fluctu- 
ations about xr in (|4I|I softens the exponent by almost 
a factor of four, as one might expect. 

The exact DOS exponent is not known for this model. 
In the next section we apply the cumulant method to the 
Tomonaga-Luttinger model, for which the exact propa- 
gator can be calculated using bosonization. 



C. Tomonaga-Luttinger model 

We consider the spinless or £7(1) Tomonaga-Luttinger 
model. The noninteracting spectrum is 



Eft = M + v F (±k - fc F ), 



(43) 



where the upper sign refers to the right branch and the 
lower to the left one. The interaction is 

V = - I dx 6ni(x)Uij 5rij(x), (44) 



8n±(x) = lira :ip±(x + a)ip±(x) :, 



(45) 



where the normal ordering is with respect to the nonin- 
teracting ground state. The matrix U has the form 



U 



We want to calculate 



Ui U 2 
U 2 Ui 



(46) 



G± (if 7i , XiTi) = -Af {T1> ± (a* , Tf )i> ± ( Xi , 7i ) e" / drV (t 5 } . 

(47) 

Make a Hubbard-Stratonovich transformation of the 
form 



e 2 



/ Sni s nj _ $D<t>-D<t> + ^gV^ 



/ D(j)_D(j) + e -H^<V<^ 



which leads to 



(48) 



J Defy- D<j> + er* 



hi four 1 



(49) 



where 
g ± (x f T f ,x i T i \(f) + ,<j>_) = 

- (T^ ± {x i T i )^ ± {x l T l )e^o dT ^ dx ^ {x ^ ]Sn ^ x ^ ) ) . (50) 
The correlation function H50|) can also be written as 

g±{x i T t ,x- 1 T i \(f) + ,(j)^) = g±(x { T{,Xin\<j>±) ■ Z T [cj) T }, (51) 
where 

and 

g±(XiTf,XiT\\<f>±) = 

- (Ttp± (xfT { )tp± {xm) e l IS dT I dx *± ix ' T) 5n± {x - T) ) 0± - 

'(53) 

Next we cumulant expand both (|52[1 and (|53|l to second 
order. For l|52|) 



Z±\<t>±] « e 



^ f drdr' f dxdx' Ii± (x — x',t— t 1 )<f>±(j)'± 



(54) 



Where fl-t is the noninteracting density-density correla- 
tion function 



U ± (x,T) = -(TSn ± (x,T) 6n±(0)) 



(55) 
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which can be written as 

n±(x,r) =G f o,±(x,r) G ,±(-a;, -r). 
The noninteracting chiral propagator is 



1 „±ikpx 

G ,±(x,t)=±— 



For (53 



2 71"? cc ± ivpT 

g±(xm,xm\<j>±) ~ Go,±(x{Tf, %in\4>±) 

^ gj dxdr Ci : ±cj}±+f dxdx drdr C2,±(p±4>'± 



(56) 



(57) 



where 



n t (fc,w) = n + (fc,w) + n_(/c,w) 



ir (lu + ik)(uj — ik) 
(65) 

The effective interaction for right movers is (with up = 1) 



?7eff(fc,u;) = f/ 4 



(cij + ik)(uj — iuk) 



where 



(58) 



where 



u= 1 + 



(w + ivk)(uj — ivk) ' 
(1/ 4 /2tt) 2 - (LW2tt) 2 

(C/4/27T) 



(66) 



(67) 



and 



1 



2tt 



(68) 



The action can be written as 
1 f dk du> 



S : 



P±(-k, -oj) U oS (k,uj) p±(k,uj). (69) 



Ci,±(a;,r) = -« — ! '- (59) 

&0 } ±{Xf,X h To) 

and Ci.± for this model reduces to 

C% t ±(x — x' , t — t') — — H±(x — x , t — t'). (60) 
Now we solve for U c e, defined by 

J dx"dr"V^{x - x", t - r")U cff (j;" - x' , t" - t') = 
5(x-x')6(t-t')1, (61) 
where 

Uf*(x,T) = 

[Ut. x {x,t) -U ± (x,T)S i+ S j+ -n^fartfi-SjJ]- 1 

(62) 

We now specialize to the DOS case where x; = Xf = 
and 1 is a 2 x 2 identity matrix. To achieve this we and assuming right movers we set p + = p. The tunneling 
Fourier transform l|61|l . This reduces (|61|) to a matrix charge density is 
equation which gives 



2 J 2tt 2tt 
The chiral tunneling charge density is 

G ,±(x { ,x,T -t) G ,±(x,Xi,T) 



p±(x,r) 



G ,±(x t ,xuT ) 



(70) 



V eS (k ) w) = ( g g 

The H — h or component is 

U eS (k, u) 



n±(fc,w) o 
o n T (fc,w) 



p(x,r) 



V F T 



27r (a; + ivFr)[x + wf{t — To) 



(71) 



(63) which satisfies 



U 4 -(Ui-U 2 2 )Tl T (k,Lu) 



(ice p(x, r) = 0(t) 0(tq — t). 



(72) 



(64) 



Fourier transforming, we find that 



p(k,uj) = — -\(e luJT ° - 1) Q(k) + (e luJT ° - e» FkT °)Q(-k)] 

iiv — vfk 



and 



p{k,u) p(-k,-u) = 



1 



(lo + ik) 2 



1 - e ^ e -kr 0+e -kr _ & -iu>T j + A . 



e(-fc) 



(73) 



(74) 
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The action therefore is 

1 f dk dio (lo + ik)(uj — iuk) U4 



S(r ) 



2 J 2-7T 2n (lo + ivk)(uj — ivk) (lo + ik) 2 

[ 1 - e^Og-fcro + g-fcro _ e -iuT \ qq^ _|_ ( 1 _ e ^o + p kr Q _ e -iWT 0( ,kT j 0( _ fc -) 



The action maybe written as S — S> + S< where 
S> 



OO oo 

— f dk [ dtu - 

87r 2 J J (to + ivk){uo — ivk)(oo + ik) 

0+ -oo 



1 — e iuJT ° e~ k ' T ° + eT kra — e~ iuT ° 



and 



(75) 



(76) 



S< = 



I dk ( duo ,^- %uk \, 

87H J J (lo + ivk)(u> — ivk)(uj + ik) 



— oo — oo 



(77) 



S> = S< under change of coordinates k — > —k and 
to — ► — w, so 5 = 2S' > . To proceed we need the large- 
tq asymptotic result 



00 p -kTo 

dk - — 

o+ k 



In r 



(78) 



where the additive constant, not shown explicitly, is cut- 
off dependent. These lead, in the large tq limit, to 



S 



u + v 



Ui r 2(1 +u) 

4tt [(1 + v)(l - v) v(l - v) 



ln(r ) (79) 



with 



\ 



Ui 

2-TTV-p 



Us 

2-KVp 



u 2 



(86) 



2izv-p 2-itvf 



Why does the cumulant method give the exact result 
for this model? The answer is that a second order cu- 
mulant expansion of the form used here is exact for free 
bosons, which are the exact eigenstates of the Tomonaga- 
Luttinger modeli^i 



S 



Ui 

4-7T 



v(l + v)_ 



Info)- 



Finally, we obtain 

S = <5 In t + const + 0(1 /to) 



and 



where 



S = 



N(e) = const x e s , 

U4 v — u (u — v)(l — v) 
4-7T v(l + v) ~ 2v(l + u) 

Ui 



1 



\ T 2irv F J \ 2ttv f ) 



2\ 1 



■ I V 27TO F / 



(80) 

(81) 
(82) 

(83) 
(84) 



2-nv-p ) V 2ttvf > 

This is in exact agreement with the bosonization result 



5 = ^-1 



(85) 



V. DISCUSSION 



Our principal result (|15|) suggests that the dominant 
effect of interaction on the low-energy DOS in a variety of 
strongly correlated electron systems is to add a semiclas- 
sical time-dependent charging energy contribution to the 
total potential barrier seen by a tunneling electron, as in 
Ref. [lj. The energy is computed according to classical 
electrostatics with a dynamically screened two-particle 
interaction. In 2D and 3D the added charge is accommo- 
dated efficiently and reaches a zero-action state at long 
times. In ID the added charge leads to diverging action, 
and hence suppressed tunneling. 

The robustness of our method has not been fully ex- 
plored, although it is known to fail qualitatively in sys- 
tems with ground state degeneracy, such as in the quan- 
tum Hall fluid. We believe the cause of this failure to 
be the non-satisfiability of the sum rule (|25|l in such sit- 
uations. In the future we plan to apply this method to 
other exactly solvable systems, such as the ID Hubbard 
and Calogero-Sutherland models. 
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